JOURNAL OF RESEARCH of the National Bureau of Standards— D. Radio Propagation 
Vol. 66D, No. 4, July- August 1962 

Magnetotelluric Fields in the Frequency Range 0.03 to 7 
Cycles Per Kilosecond: Part I Power Spectra 

C. W. Horton and A. A. J. Hoffman 

Contribution from Department of Physics, The University of Texas and Defense Research 

Laboratory, The University of Texas 

(Received November 29, 1961; revised February 1, 1962) 

Power spectra of the horizontal components of the magnetic field and the telluric field 
are computed for data recorded on 1 and 2 September 1957 at the Soviet Magnetic Observa- 
tory in Tbilisi. Power spectra of the East- West telluric field component are computed for 
20 September 1957 for Soviet stations located at Lvov, Tbilisi, and Ashkhabad. All analyses 
are based on microfilm copies furnished by the IGY World Data Center A. Each power 
spectrum shows a frequency dependence over a frequency range of 0.03 to 7 cycles per 
kilosecond of the form Po/~" where n varies from 1.1 to 2.5. The coherency between or- 
thogonal components of the telluric and magnetic fields is computed for the data from Tbilisi. 
Magnetotelluric power spectra from USSR, Canada, Texas, and Massachusetts are plotted 
on a common graph to show the frequency dependence over the range 0.03 to 400 cycles per 
kilosecond. 



1. Intrcxluction 

The aim of the present study is the computation of 
power spectra of the magnetotelluric fields at the 
earth's surface. One of the authors ((/WH) spent 
a week in August 19G0, at the IGY World Data Cen- 
ter A, Washington, D.(\, examining the microfilm 
copies of data available at that time. The available 
data were examined with a view to finding simul- 
taneous recordings of telluric and magnetic fields at 
one station so that Cagniard's theory of magnetotel- 
hu^ic prospecting [Cagniard, 1953] could be applied 
to the earth. A second goal of the study is the com- 
parison of simultaneous recorcHngs at different sta- 
tions. With these aims in mind tlie microfilms were 



searched for days on which the fluctuations were 
large so that the records could be read easily but 
which were free of storms and sudden commence- 
ments so that the records would approximate sta- 
tionary time series. This last condition was impor- 
tant since the data would be subjected to spectral 
analysis. 

Of the microfilms available at the time of the visit 
to the Data Center those for the month of September 
1957, for the IGY stations Lvov, Tbilisi, and 
Ashkhabad,^ USSR, were most suitable for the 
present purpose. Figure 1 shows the location of the 
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Figure 1. The location of the three Soviet IGY Stations Lvov, Ashkahbad, and Tbilisi. 
Base map copyright by Univ. of Chicago. 
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stations, and figures 2 and 3 sliow typical samples of 
the data. The reader is referred to the publications 
of the I.A.G.A. [1957] for detailed descriptions of 
these observatories. 

The magnetic and telluric fields recorded at Tbilisi, 
on 1 and 2 September 1957, were subjected to fairly 
extensive analysis. The East- West component of the 
telluric fields at the three stations were analyzed for 
simultaneous recordings made on 20 September. 
This gave a cross section in space on one day and a 
sampling in time at one station for one of the com- 
ponents of the telluric field. 

The telluric recordings give the East-West and the 
North-South components directly. The magnetic 
recordings, however, give the horizontal intensity H 
and the declination D. These values were converted 
to North-South and East-West components for the 
present analysis. 

2. Methods of Processing the Data 

After suitable days were selected, photographic 
enlargements were made of the records and the ampli- 



tudes of the traces were read at equally spaced times. 
The time interval, A^, between successive samples was 
chosen small enough to avoid loss of the higher fre- 
quency components. The data were processed by the 
methods described by Blackman and Tukey [1959] 
whose notation and terminology will be followed 
closely. Their procedures have been used by other 
workers in this field [Munk, Snodgrass, and Tucker, 
1959; Cantwell, I960]. Table 1 contains a summary 
of the parameters associated with each sample of data 
discussed in this paper. 

After the ordinates were read, they were subjected 
to an arithmetic operator [Horton, 1962] which acted 
as a high-pass filter. The power spectrum of the 
filtered values is more nearly constant than that of 
the unfiltered values so the filter serves as a ^^pre- 
whitening" filter. The N filtered values yi were 
used to compute the mean lagged products 



Cr=j^^T.^yiyi^r,r=0,l,2, . . .,30. (1) 
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Figure 2. A typical sample of the luagnetograms analyzed. 
Record for Tbilisi 2 September 1957. 














Figure 3. A typical sample of the tellurograms analyzed. 
Record for Lvov 20 September 1957, 
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Table 1. Summary of the data presented in this paper 



Sample No 

Date 

Telluric field 

Magnetic field 

Paper speed 

Length of sample hr__ 

At, sec 

No. of samples 

Decimation 

New At sec__ 

No. of decimated samples 

Plotted in figure 



20 Sept. 57 
E-W 

X 

22 mm/hr 
24 
150 
579 
4 
600 
144 
5 



20 Sept. 57 
E-W 

X 

90 mm/hr 

5.5 

31.6 

624 

4 

126.4 

156 



20 Sept. 57 
E-W 

X 

78 mm/hr 

7.8 

47.3 

590 

4 

189.2 

147 



2 Sept. 57 

E-W 

N-S 

22 mm/hr 

10 

72 

501 

4 

288 

125 

8.9 



2 Sept. 57 

N-S 

E-W 

22 mm/hr 

10 

72 

501 

4 

288 

125 

10.11 



31 July 59 
E-W 

N-S 

12 in./miu 

1.25 

0.75 

3, 030 

4 

3 

757 

13,14 



The coefficients Cr represent approximations to the 
autoco variance function whose Fourier (cosine) 
transform yields the power spectrum. Thirty-one 
estimates of the power spectral density are obtained 
from the following cosine transform 



Vr = M [C^o 



+2 1] C, cos {(irT/30) + Cso cos rr 



} 



(2) 



The resulting power spectral estimates are averaged 
in groups of three to give improved values. This is 
the process that Blackman and Tukey refer to as 
hanning. Finidly, the effects introduced by the 
high-pass filter nre removed to obtain the values 
plotted. 

In order to increase the resolution at lower 
frequencies, one must either increase the range of 
the index r of eq (1) or increase the sampling interval 
A^. Before one can safely increase A^, a low-pass 
filter must be used to insure that no significant 
power occurs at frequencies greater than the Nyquist 
frequency 1/2 A^. When the power spectrum is 
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Figure 4. Frequency response {amplitude) of the low- and high- 
pass filters used in processing the data. 



contaminated by power above the Nyquist frequency 
one sometimes describes this effect as ^ ^aliasing.'' 
The second alternative mentioned above, which is 
preferable in that it requires less computation, was 
followed. The sampled data were filtered with a 
low-pass numerical operator and decimated by 
selecting only every fourth value. After the power 
spectral density estimates were computed, the values 
were corrected for the effect of the low-pass filter. 

The frequency responses (amplitude) of the two 
filters are shown in figure 4. Although the abscissa 
of the curve is labeled for a sampling interval of 
A^ = 72 sec, it can be corrected to any other sampling 
interval by a simple change of scale. The low-pass 
filter is a symmetric 9-terin operator given by Munk, 
Snodgrass, and Tucker [1959]. The high-pass filter 
is an unsymmetric 18- term operator designed by 
Horton [1962]. Since the low-pass operator is sym- 
metric, it does not produce a phase distortion. The 
phase delay of the high-pass operator has no effect 
on the powder spectra. As one may see from figures 
2 and 3 tlie data do not liave any d-c drifts and no 
correction for such drift was necessary. 

3. East- West Telluric Components for 
20 September 1957 

Figures 5, 6, and 7 are gi'aphs of the power sptH'tra 
of the East- West component of the telluric field at 
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Figure 5. The power spectrum of the East-West telluric field at 

Ashkahbad 20 September 1957, 

80 percent confidence is indicated. 
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Figure 6. The power spectrum of the East-West telluric field at 

Lvov 20 September 1957. 

80 percent confidence is indicated. 
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Figure 7. The power spectrum of the East-West telluric field at 
Tbilisi 20 September 1957 

80 percent confidence is indicated. 

Ashkhabad, Lvov, and Tbilisi, respectively, for 20 
September 1957. The graphs are self-explanatory 
and require very little comment. Each of the power 
spectra obtained from the low-pass filter shows some 
evidence of aliasing. 

The three power spectra can be approximated 
with a formula 



p{i)=pj- 



(3) 



with "a value of n between 1.5 and 1.6. Actu- 
ally this formula is not a satisfactory fit since 
one would select n nearer one for the lower decade 
of frequency and nearer two at the upper decade of 
the spectrum. The question of an empirical formula 
will be considered further in section 6. 



It is interesting to compare the power density at 
each of the stations at a frequency of 1 c/ks. If one 
uses Lvov as the reference, the power densities at 
Ashkhabad and Tbihsi are —6.3 db and —18.8 db, 
respectively. These relative values undoubtedly 
reflect the underlying continental structure. For 
example, the topography at Tbilisi is higher and more 
rugged than the other two stations. The moderately 
saline Black Sea (salinity 17.6 to 22.2 %o) is 
between Lvov and Tbilisi but nearer the latter. 

4. Magnetotelluric Spectra at Tbilisi 

As mentioned above the records made at Tbilisi 
on 1 and 2 September 1957 w^ere subjected to a 
detailed analysis. The reader is referred to the 
Ph. D. thesis of Hoffman [1962] for more graphs and 
tabular values of the power spectra Most samples 
were 10 hr long. A few computations w^ere made for 
20-hr samples but these did not show enough 
difference from the 10-hr samples to warrant their 
presentation. For some reason that cannot be 
discovered, when the East-West component of the 
magnetic field and the North-South component of 
the telluric field were analyzed with the aid of the 
low-pass filter, the power spectra were erratic and 
unreliable. This showed up both in numerous 
negative values for the power spectral densities but 
also in ridiculously large values in the coherency. 
These low-passed data w^ere rejected and the re- 
sponses of the high-pass filters only are used in this 
paper. 

The individual magnetic and telluric power spectra 
for 2 September are illustrated in figures 8 through 
11. The spectra of the North-South magnetic and 
the East- West telluric components do not show any 
striking departures in behavior from the curves for 
20 September. The empirical formula in eq (3) can 
again be made to fit the magnetic spectra with n=2 
and the telluric spectra with n between 1.1 and 1.2. 
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Figure 8. The power spectrum of the North-South magnetic 

field at Tbilisi 2 September 1957, 

80 percent confidence is indicated. 
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Figure 9. The power spectrum of the East-West telluric field at 

Tbilisi 2 September 1957. 

80 percent conftdonce is indicated. 
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Figure 10. The power spectrum of the East-West magnetic 

field at Tbilisi 2 September 1957. 

80 percent confidence is indicated. 
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Figure 11. The power spectrum, of the North-South telluric 

field at Tbilisi 2 September 1957. 

80 percent confidence is indicated. 



If we compare power densities at 1 c/ks and use 1 
September as the reference, tlie power levels of the 
East-West telluric field are db, +0.6 db, and —0.4 
db on 1, 2, and 20 September, respectively. The 
power levels of the North-South magnetic fields are 
db and +0.6 db on 1 and 2 September, respcMMivf^ly. 

5. Coherency of the Magnetotelluric Field 
at Tbilisi 

It has been known for some time [Chapman and 
Whitehead, 1923] that any horizontal component 
of the telluric field is related causally to the orthogonal 
horizontal component of the magnetic field. The 
degree of correlation between two such components 
can be evaluated from the coherency defined as 



Coh^^(f) = 



\^eh{^)\ 



^^EE{mHH{i) 



(4) 



where ^eh(J) is the cross spectrum between the two 
field components and ^^^^(f) and ^hh{^) are the 
corresponding powei- spectra of the individual 
series. The cross spectrum is approximated from the 
saniphnl data with the aid of the lagged cross products 



1 N-r 

iV — V j=o 



r-0,±l, 



±m. (5) 



When X and y are dift'erent time series, C-r^C+r, 
in general, and the Fourier transform is complex. 
Therefore eq (2) must be supplemented by a similar 
expression with sines. Thus, the corresponding 
power spectral density estimates ^eh ^i'^' complex 
and one must use absolute values in eq (4). Figure 
12 contains a few plots of the coherency between 
components of the magnetotelluric fields at Tbilisi 
on 1 and 2 September. 




E-W TELLURIC AND N-S MAGNETIC 
I SEPT 1957 




E-W TELLURIC AND N-S MAGNETIC 
2 SEPT 1957 
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N-S TELLURIC AND E-W MAGNETIC 
I SEPT 1957 



Figure 12. The coherency function CohEnif) for the magneto- 
telluric field at Tbilisi 1 and 2 September 1957. 

The Nyquist frequency of each spectrum is 6.94 c/ks except for the upper two 
curves on the right. The Nyquist frequency for these two curves is 1.74 c/ks. 
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6. Comparison With Other Magnetotelluric 
Data 

Althougli each of the spectra plotted m figures 5 
through 11 cover 2.4 decades of frequency, it is evi- 
dent that this is a small section of the total mag- 
netotelluric spectrum. In order to obtain a more 
comprehensive picture of the magnetotelluric spec- 
trum selections from the present data and other 
spectra available to the authors are combined on a 
composite plot in figures 13 and 14. 

Since the Canadian data presented in these fig- 
ui-es are new data that have not been published 
before, a brief description of the experimental ar- 
rangement is given. In the summer of 1959 Sir 
Charles Wright, Dr. J. H. Duffus, and other staff 
members of the Pacific Naval Laboratory, Esqui- 
mault, Canada, established a field station at the 
Suffield Experimental Station near Medicine Hat, 
Alberta, Canada, to measure magnetic field fluctua- 
tions. Dr. K. Vozoff of the Universit3' of Alberta 
set up earth current electrodes at this station for 
the measurements of telluric fields. Mr. Jerry 
Joyner, an employee of Texas Instruments, Inc., 
Dallas, was a member of the field party. Mr. 
Joyner very kindly supplied the authors with simul- 
taneous recording of the East-West telluric field and 
of the North-South magnetic field. The component 
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A TBILISI, USSR I SEPT 1957 
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Figure 13. A composite graph cf the power spectra of the 
North-South component of the magnetic field for USSR, 
Canada, Massachusetts and Texas. 

An absolute calibration of the Canadian data is not available so these points 
can be displaced vertically an arbitrary amount. 



of the magnetic field was measured with a metal- 
cored solenoid so that the recording showed the 
time derivative of the magnetic field. The com- 
puted power spectrum was corrected to elimmate 
the eft'ect of this derivative. The telluric recording 
showed the potential difference between two copper 
electrodes separated 1,000 ft. The paper speed of 
the recorder was 12 in. per minute and the records 
were sampled at an interval A/=0.75 sec. A sam- 
ple of the data was analvzed which started 0424 Z, 
31 July 1959 Greenwich' time or 9:24 p.m. 30 July 
local time and lasted 75 min. The results of the 
low-pass filter only are included in figures 13 and 
14. No absolute sensitivities are available to the 
authors so the vertical positions of these two curves 
were adjusted to fit the other data; the frequency 
dependence is still meaningful, however. 

The magnetotelluric data for Littleton, Mass., 
were recorded by Dr. Thomas Cantwell on 5 Decem- 
ber 1959. These spectra are given in his thesis 
[Cantwell, 1960] as Case 27-3. 

The magnetotelluric data for Austin, Tex., were 
recorded 27 March 1961 by Dr. H. W. Smith and 
Mr. F. X. Bostick of the Electrical Engineering 
Research Laboratories, The University of Texas. 

The locations of these various stations are shown 
in figure 15. Despite the fact that the data were 
recorded over a period of 2^ years and extend over 
a longitude range of 170°, the data appear to fit a 
uniform trend. It is suggested that the trends es- 
tablished in figures 13 and 14 are a first approxima- 
tion to the ambient magnetotelluric field at the 
earth^s surface. Although fluctuations of the power 
densit}' at one station may vary 10 db (Dr. H. W. 
Smith, personal communication) and the level be- 
tween two stations may be as great as 20 db (Lvov 
and Tbilisi), these variations are small compared to 
the variation of 80 db in the telluric field plotted 
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Figure 14. A composite graph of the power spectra of the 
East-West component of the telluric field for USSR, Canada, 
Massachusetts, and Texas. 

An absolute calibration of the Canadian data is not available so these points 
can be displaced vertically an arbitrary amount. 
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FiGUHE 15. A world Diap showing the locations of stations referred to in fujures IS and 14- 
Base map copyright by Univ. of Chicago. 



in fig:ure 14. The corresponding" variation in the 
power density ot* the magnetic field plotted in figure 
13 is 150 dh.* 

The following empirical formulas represent useful 
approximations to the data presented in figures 13 
and 14. When / is expressed in cycles per kilo- 
second, the magnetic field can he expressed hy 

(gammas)Vc/s, 6.1 c/ks</<400 c/ks, ... (6) 

Superimposed on this background spectrum is an 
intense band of power at 80 c/ks (the frequency of 
pc's) and a smaller one at 180 c/ks. There is slight 
evidence of the power baTid at 80 c/ks in the (Cana- 
dian data. The electric field can be expressed by 
(again / in c/ks) 

P^(f)-=10^[l + 100/+20y2 + f]-i (mv/km)Vc/s, 

0.03 c/ks</<400 c/ks, ... (7) 

The experimental data for P^(f) scatter more strongly 
* than the data for Ph{^)- This is undoubtedly due 
to the influence of the earth's structure on the 
telluric field. 
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